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Abstract 

The purpose of this article is to discuss cluster expansions in dense 
quantum systems as well as their interconnection with exchange cycles. 
We show in general how the Ursell operators of order I > 3 contribute 
to an exponential which corresponds to a mean-field energy involving the 
second operator U2, instead of the potential itself as usual - in other words, 
the mean-field correction is expressed in terms of a modification of a local 
Boltzmann equilibrium. In a first part, we consider classical statistical me- 
chanics and recall the relation between the reducible part of the classical 
cluster integrals and the mean-field; we introduce an alternative method 
to obtain the linear density contribution to the mean-field, which is based 
on the notion of tree-diagrams and provides a preview of the subsequent 
quantum calculations. We then proceed to study quantum particles with 
Boltzmann statistics (distinguishable particles) and show that each Ursell 
operator U n with n > 3 contains a "tree-reducible part", which groups 
naturally with XJ% through a linear chain of binary interactions; this part 
contributes to the associated mean-field experienced by particles in the 
fluid. The irreducible part, on the other hand, corresponds to the effects 
associated with three (or more) particles interacting all together at the 
same time. We then show that the same algebra holds in the case of Fermi 
or Bose particles, and discuss physically the role of the exchange cycles, 
combined with interactions. Bose condensed systems are not considered 
at this stage. The similarities and differences between Boltzmann and 
quantum statistics are illustrated by this approach, in contrast with field 
theoretical or Green's functions methods, which do not allow a separate 
study of the role of quantum statistics and dynamics. 
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1 Introduction 



A widely used formalism in quantum statistical physics is the formalism of 
Green's functions H, ^j, where the techniques of second quantization and 
field theory are used from the beginning; the notion of exchange operators of 
indistinguishable particles is of course contained, but in a completely implicit 
way. In the formalism of Ursell operators ^, M, [|], the starting point is first 
quantization with numbered particles, so that the role of exchange cycles be- 
comes completely explicit: these cycles appear clearly in all diagrams and, for 
instance, they are the only source of diagrams for the ideal quantum gas. This 
reduces the distance between the formalism and, for instance, numerical calcu- 
lations such as the PIMC method (Path Integral Monte Carlo), where particles 
are also numbered and the exchange cycles are explicitly sampled by random 
choices. Moreover, it becomes possible to assume that the particles obey Boltz- 
mann statistics, just by "switching off the cycles" , a task that would be difficult 
in the Green's function formalism. Needless to say, this does not mean that 
Green's functions are, in general, less powerful than the Ursell formalism! The 
opposite is actually closer to reality: for instance, Green's functions handle time- 
dependent problems easily, while this is not the case in the Ursell formalism. 
But it remains true that, if one is interested in a detailed discussion of the effects 
of quantum statistics, it becomes more straightforward to resort to the Ursell 
formalism. 

In this article, we consider dense systems, for which it is not necessarily 
possible to limit oneself to first order density effects. In contrast to the situation 
in a dilute gas, a given particle may interact frequently with several others at 
the same time, and even liquefaction may take place. Therefore, we will no 
longer ignore all Ursell operators beyond U2, as was done in most of previous 
work in this formalism; operators C/3, U4, etc. now become important. One 
may actually wonder what the role of these higher order operators is in general, 
and why exactly it is possible to ignore their role in a dilute system, as was 
done in J6| for instance. We will see that part of their contribution (what we 
will call their tree reducible part) groups naturally with U2 through a chain of 
binary interactions, and builds an exponential of the mean-field energy. In other 
words, instead of making the problem more complicated, this contribution of 
the higher order operators builds exactly the exponential of U2 that is needed 
to reconstruct a simple and natural expression of the mean-field. Nevertheless, 
this mean-field is expressed in terms of the matrix elements of U2, instead of 
the usual matrix elements of the potential itself; in a sense, what we obtain is 
the exponential of an exponential, since U2 itself contains exponentials of the 
Hamiltonian and corresponds physically to the local change of the Boltzmann 
equilibrium. As a consequence, the logarithms that appeared in ||, and had to 
be expanded to first order in density, are actually spurious - in other words, 
this first order expansion was actually unnecessary. In addition, the rest of the 
contribution of the higher order Ursell operators, the irreducible part, vanishes 
unless three particles (or more) are all close together, and is really characteristic 
of many-body collisions and of dense systems. 



2 



2 Classical statistical physics 



In this part, we quickly review the classical cluster expansion for the parametric 
equation of state of a fluid; we first recall the results of classical statistical 
mechanics as a point of comparison. 



2.1 General formalism 

The general expression of the equation of state was derived by Mayer [Q and 
Ursell see also the books by Uhlenbeck and Ford || and Hansen and Mc- 
Donald [1CJ. A classical system of massive particles, with mass m, is supposed 
to be contained in a box of finite volume V, with periodic boundary conditions 
(translationally invariant system); their Hamiltonian is the sum of the kinetic 
energies plus the interaction energy, which is the sum over all pairs of particles 
of the interparticle pair potential Vij = V(ri — rj). An useful function in the 
calculation is the Mayer function defined by: 

/y = exp(-pVij) ~ 1 (1) 

This function is everywhere bounded and goes to almost zero when the distance 
between particles i and j is much larger than the range of the pair potential; it 
is the classical equivalent of a second Ursell operator U%(i,j) ||. 

The classical cluster expansion is the expansion, in the grand canonical en- 
semble, of the pressure p and of the density p of the gas in series of the fugacity 
z = e***, where /j, is the chemical potential and (3 = l//cgT is the inverse 
temperature. If A is the de Broglie thermal wavelength, the expansion can be 
written as: 

with: 

(N) d f\nZ g 



1=1 



where Z gc (f3, V,z) is the grand canonical partition function; the definition of 
the cluster integrals bi(T, V) is the same as that in the book of Mayer & Mayer 
[[7j. The first coefficients are given by: 

h(T,V) = -^Jd 3 n = l (4) 



2!V 



d 6 n d A r 2 f 12 = - j d A r 12 /(r 12 ) 



(5) 



h{T, V) = -^J d 3 n d 3 r 2 d 3 r 3 [f 12 .f 13 + f 12 f 23 + f 13 f 23 + Aa/ia/as] (6) 
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and so on for / > 4; these definitions include a 1/V factor so that bi(T, V) 
remains intensive when the volume tends to infinity. 

It is useful to introduce diagrams to represent the integrals. A /-particle 
cluster diagram is made of I numbered circles, representing the particles, be- 
tween which one draws any number of lines (also called links, and representing 
the fi/s), each line joining distinct pairs of circles. The diagram is said to be 
connected when one can go from any particle to any other particle in the cluster 
following the lines. The general definition of the cluster integrals is then: 

bi(T,V) = yp^j x sum over all possible Z-particle clusters (7) 
For example: 

6 1 (T,V) = -i-(.) = l (8) 



M^l^H (9) 

and: 

h(T,v) = ^ (A + A + A + A) (io) 

2.2 Irreducible diagrams, energy shift and equation of 
state 

It turns out that, defined in this way, diagrams carry redundant information: 
as soon as / > 3, most of the diagrams contributing to bi can be constructed 
as products of smaller diagrams (already contained in 6 ; < with I' < I). For 
instance, this is the case of all diagrams in the right hand side of except 
the last. It then becomes convenient to define the notion of irrcducibility: an 
irreducible cluster is such that, if any line (which is essentially a 62) is removed 
from it, it never splits into disconnected clusters; conversely, in a reducible 
cluster, it is possible to find at least one line that, when cut, will decompose the 
result into two separate clusters (since each line is essentially a 62, reducibility 
and irrcducibility are defined here in terms of binary coefficients 62 's, but more 
general definitions are conceivable). 

Irreducible diagrams play a special role if one introduces an exponentiation 
of the fugacity expansion of the density @). Let us define an energy shift A, 
which shifts all energy levels e&, by the relation: 

p = e-« A -">i (11) 
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It can then be shown (see |yj, @, @, Q) that A is given by the series: 

oo 

-/3A = ^>(T,V)p' (12) 

where the new coefficients (3i are defined by: 

Pl-i(T, V) = jz -rr— x sum over all irreducible Z-particle clusters (13) 

or, equivalently: 

A_i(T, V) = I X Irreducible part of b t (T, V) (14) 

For example: 

Pi(T, V) - i (•— ■) = J d 3 r 12 /(n 2 ) (15) 

and: 

A(T ' V) = 2Fv(A) 

One can also show that irreducible clusters are directly related to the equation 
of state (at finite volume), from which the variable z has been eliminated: 

f "i-Eyxi^v)^ (it) 
9 i=x ti " x 

The proof of this result is given in the book by Mayer & Mayer H. Another 
method of calculation was introduced by Van Kampen jL5|, who works with 
the canonical ensemble instead of grand canonical; see also |l(|. In these cal- 
culations, the thermodynamical limit is taken before the I summation while, 
here, we do not take this limit; the reason is that, as discussed by Lee and 
Yang [0, taking directly the limit for each coefficient reduces the validity of 
the calculation to the gaseous phase only. 

2.3 Tree-reducible diagrams; mean-field 

Let us now suppose that, among all clusters contributing to 6;, we keep only 
those having I particles and / — 1 links. They correspond to the minimally con- 
nected diagrams: if any link is removed, these diagrams split into two different 
diagrams. They are, not only obviously reducible (as soon as I > 3), but also 
"fully reducible" , since the removal of any line splits the diagram into discon- 
nected clusters. For instance, the fully reducible part of is: 

b§(T) = j d 3 n d 3 r 2 d 3 r 3 [f 12 f 13 + fa fa + fa fa] (18) 
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Figure 1: A tree-diagram corresponding to a 10-particle cluster; the tree sym- 
bolizes one of the terms that are contained in the tree-reducible part of 610. The 
construction method that is used to build the branches from the structure of 
the integral is explained in the text. 



from which the last term that appear on the right hand side of (|6|) has been 
eliminated. Graphically, we will represent these clusters by diagrams that have 
the structure of a tree, such as that of figure [j], and for this reason we shall call 
them "tree-reducible diagrams" . 

In order to construct this tree-diagram, we have to choose one numbered 
particle as the "root particle" , for instance particle 1 . Then, among all particles 
j that are connected to particle 1 by a link f±j , we select that with the smallest 
value of j; let us call k this particle. The same operation is then made again 
from particle k: one identifies the particle m to which it is related which has 
the smallest numbering, and puts this particle as the next in the branch; the 
same operation goes on until one reaches the end of the first branch, which is 
drawn vertically by convention. One then goes back along this branch towards 
the origin and identifies the first particle which is connected to another particle 
outside of the branch; this particle is a branching point, the source of another 
branch, which is built in exactly the same way, and drawn directly on the 
right of the first branch. The same process continues and one adds successive 
branches, which are drawn in a clockwise order from their source; it stops when 
all particles are included in the tree. This construction provides a well defined 
geometry for the tree-diagram but, clearly, different tree-diagrams may have the 
same numerical valuejj]. It is nevertheless convenient for our reasoning below to 
classify their contributions according to the various structures of trees. 

For instance, the reducible part of &3 then becomes: 



1 In fact, the value of these diagrams depends only on the total number of nodes (i.e. 
particles) I (see Appendix I). 




(19) 
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(from now on, we assume that the factor 1/V is included in the value of the 
tree-diagram, in order to simplify the notation). 

Let us now consider a given tree-diagram and evaluate its weight; the ques- 
tion is to determine how many different numberings are compatible with the 
rules that we have used to build it. Suppose then that we put random numbers 
into the I — 1 nodes that are available in the tree. The result will be acceptable 
only if a correct clockwise ordering of the particle numbers is obtained at each 
bifurcation node. Let us call n, ra, .. the branching factors (or ramification 
factors) of the nodes, i.e. the number of secondary branches at each node; a 
linear diagram has only r = 1, a diagram with one binary bifurcation has one 
r = 2, and r takes the values 3, 4, etc. when more branches start from the 
same source. Each node will then introduce a probability l/r^ ! for this correct 
ordering to be obtained. Since there are (I — 1)! ways of distributing the I — 1 
particles among the nodes, the final result is that the weight of the tree-diagram 
is: 



ILr<! 



(20) 



If we call Tdiag the value of a tree-diagram (including the 1/V factor), we can 
then express the tree-reducible I particles cluster as: 



{tree diagrams} 



Tdi, 



(21) 



Inserting this expression into the z-expansion of the density ([s|) , we notice that 
various simplifications take place in the coefficients, so that we obtain: 



z 1 z/A 3 T (z/A 3 ) 
A 3 | lF X ^2P 



( 



2! 



V 




2! 

IT 



(22) 



It is now possible to regroup diagrams according to the value of the branching 
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factor ri at their root: 



^ 1 + Tl 



+ (£) 



(23) 



One then notices a sort of self-similarity property of the expansion: at each 
secondary nodes, one gets an expansion that provides again all the terms that 
are contained in p. Therefore: 



P = 



2! 



(24) 



We now recognize the development of an exponential: 

P=^{l + /3iP+^x [p lP f + ... + -L x [/3 lP p + . 

where (3\ — 2&2 is defined in (|i~5|), and finally obtain the very simple result: 

p=^exp(/3 lP ) (25) 

or: 



-j9A = ftp 



(26) 



Therefore, when only tree-reducible diagrams are taken into account, the energy 
shift becomes exactly proportional to the density; we recover what is usually 
called the mean-field approximation. The method we have used is different from 
the method of refs. jll, [l2], [l3| it relics on the branching properties of tree- 
reducible diagrams and not on the use of complex variables. Moreover, as we 
will see below, the method can be transposed to quantum mechanics. 

Of course, equation ( p6| ) can also be obtained by the same method as these 
references, in a particular case: one assumes that the only non-zero irreducible 
cluster integral is (3±. Clearly, in this situation, all cluster integrals reduce to 
their tree-reducible value bf-, which can be expressed as a product of terms 62 = 
(3i/2 (sec Appendix I): 



R _ 



-Ni ft 



1-1 



(27) 



8 



where Ni is the number of clusters with I numbered particles and I — 1 links, 
i.e. the number of classical tree-diagrams containing I particles: 

{tree diagrams} 

If we insert the corresponding value of 6; into equations ^ and (||) , we obtain 
again equation (p6|). Moreover, equation then becomes: 

fo=l-\fc(T,V) p = l-b 2 (T,V) p (29) 
P * 

which contains only a linear density correction. 

3 Quantum statistical physics: Boltzmann statis- 
tics 

We now leave classical statistical mechanics and reason within quantum mechan- 
ics; in a first step, we consider distinguishable particles (Boltzmann particles), 
postponing the discussion of the effects of quantum statistics to the next sec- 
tion (bosons or fermions). We first introduce our notation and then discuss the 
introduction of exponentials from the structure of tree diagrams. 

3.1 Notation 

We assume that the Hamiltonian of the system of N non-relativistic particles 
is: 

N 

i—l i<j 

where Hq(z) is the one-particle energy (sum of its kinetic energy plus coupling 
to an external potential) and where the binary interaction potential Vij is a 
function of the distance between particles i and j: 

Vij = v (\ r i - r j\) = v ( r ij) (31) 
As in pL we define the Ursell operators: 

E/i(l)=exph9flb(l)] (32) 

and: 

U 2 (l, 2) = exp [-13 [flo(l) + #o(2) + V 12 ]} - U 1 (1)U 1 {2) (33) 

and so on for higher order Ursell operators Us(l,2, 3), (74(1,2,3,4), etc. All 
the t/;'s for I > 2 have a clustering property; for instance, the diagonal matrix 
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elements of < ri, r 2 | U 2 | ri, r 2 > tend towards zero when the distance | ri — r 2 | 
becomes larger than some microscopic distance; in other words, the matrix 
elements of the same operator in the momentum representation are proportional 
to the inverse volume 1/V. One can then express the canonical partition function 
Zn as a sum of products of traces of operators Uj 's, each corresponding to a given 
diagram. Here, since we are dealing with Boltzmann particles, the number of 
diagrams is smaller than in [Q ; exchange cycles do not occur, which is equivalent 
to limit them to cycles containing one particle only - in other words we exclude 
all those diagrams that contain horizontal lines^j. Going to the grand canonical 
ensemble, one can then show that the logarithm of the corresponding partition 
function Z gc , is given by the following sum of traces: 

00 1 

lnZ g . c . = J2^Tr 1 ,2,...,i{U l (l,2,...,l)} (34) 
1=1 

where the factor 1/Z! corresponds to the weight of the diagram, which arises 
because there are 11 equivalent ways to distribute I numbered particles inside 
the Ursell operator, as discussed in 

Rather than studying a thermodynamic potential, it is often more convenient 
to focus the discussion on the single particle density operator p/(l), as in ref. 

One then gets: 

PI (1) = zUx{l) + z 2 Tr 2 {U 2 (l, 2)} + ^ Tr 2 , 3 {U 3 (l, 2, 3)} + ... (35) 

Now, because particle 1 is "tagged" , it plays a special role, so that the generic 
term of this series is: 



Ojj% ! {^i(1.2, 3, ...,!)} 



(36) 



with a weight 1/(1 — 1)! corresponding to the equivalent distributions of all 
the untagged particles. In the usual graphical representation or Ursell diagrams 
[Q 0], where vertical lines symbolize Ui operators (for I > 2) and horizontal lines 
exchange cycles (we have already mentioned that, since here we are dealing with 
distinguishable particles, no cycle of length greater than one occurs), equation 
B~5|) becomes: 



pi(l) = 1—1 + i-IL, -(- — X h 



(37) 



At this stage, it is convenient |], |T^] to express the second rank operators U 2 in 
terms of the operator U 2 defined by: 



C/ 2 (l,2) = y/Ui (l)J7i(2) U 2 (l,2) y/ui (XjUi (2) (38) 



2 As usual, we retain the overall 1/7V! factor in the symmetrizer, in order to avoid some 
well-known difficulties of classical statistical mechanics (Gibbs paradox, etc.). 
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(for the moment, when I > 2, we do not specify the exact relation between 
an Ui operator and Ui, but we will come back to this point later). We note 
that, here, U 2 is defined symmetrically, with square roots of U\ operators on 
each side (U\ is a positive operator), so that U 2 is Hermitian - this was not 
the case in [[l9|. The introduction of U 2 has two advantages. First, factorizing 
the kinetic energy brings the formalism closer to classical statistical mechanics, 
where kinetic energy always factorizes out exactly; in other words, U 2 (i,j) is 
the quantum equivalent of /y , although it is not strictly analogous (in quantum 
mechanics operators do not necessarily commute, so that the kinetic energy 
can still play a role in 1/2(1, j))- Second, in the limit of low energies, it is 
possible to make use of the MIME (momentum independent matrix elements) 
approximation, where the diagonal matrix elements of U2 are constants - see 
for instance the characterization of the diagonal elements of U2 in terms of the 
Ursell length in pJJ. Graphically, in order to distinguish [/2's from U 2 's in 
diagrams, we use two vertical parallel lines in the former case, as in 

§ and g, 

but only a single line for [/2's. 

3.2 Tree-reducible part of U\ and V\ 

For any value of I, we now define the tree-reducible operator (or fully reducible 
2^ 

operator) l/ x by analogy with the tree-reducible part of a classical cluster 

Since the classical links correspond to U2 operators in quantum mechanics, and 

since the minimal number of links is I — 1, we will express the fully reducible 

operator Ui as a product of I — 1 operators U 2 ', in addition, we also have to 

take into account the fact that operators do not necessarily commute with each 

j£ 

other and apply an appropriate symmetrization. We will therefore define V \ by 
a double sum: 

Ff(l, 2, ...,/)= ^-y £ U 2 (.,.)U 2 (.,.)...U 2 (.,.) (39) 

{°P- orderings} 

The first sum symbolizes all different ways to choose sets of / — 1 operators U 2 so 
that the product correspond to a minimally connected cluster - in other words, 
in all classical tree-diagrams, we replace each /y by the corresponding 1/2(1, j)- 
The second sum, together with the factor 1 / (/ — 1 ) ! , corresponds to an average of 
the product of all operators over all possible orderings of these / — 1 operators. 
With this definition, U l is obviously Hermitian as well as symmetrical with 
respect to all particles 1, 2,..., I. For instance: 

U*(l, 2, 3) = ±[U_2(1, 2)^2(1, 3) + U_2(l, 3)^2(1, 2)] 

+ ^2(1,2)02(2,3) + U 2 (2,3)U 2 (1,2)} (40) 
+^2(1, 3)C/ 2 (2, 3) + U 2 (2, 3)U 2 (1, 3)] 

We now wish to define an operator Uf 1 that we will call the tree-reducible 
part of Ui . The first idea that comes to mind is to mimic ( |38|) and to define Uf 1 
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Figure 2: In this diagram, the single lines symbolize XJ i operators (as opposed 
to double lines which symbolize ?7 2 's in Ursell diagrams). The figure illustrates 
symbolically how y/Ui's (symbolised as squareroots in the left part, as rect- 
angles in the right part) are inserted into each product of J7 2 's that appears in 
U l , in order to build an operator 



by just multiplying E/j on both sides by a product of yUi's. This is possible, but 
it turns out that this definition would be less convenient than another slightly 
different possibility, where some square root operators are inserted at different 
places. To introduce this definition, we first remark that each term in the first 
summation of ( |39| ) can be associated with a tree, exactly as in the classical tree 
of figure [I], where particle 1 is put at the root, and all the rest of the diagram 
is built exactly in the same way (the /y's are replaced by /7 2 (i,j)'s). The 
problem now is how to properly "dress" this diagram with y/Ui's in order to 
build a proper Uf 1 ; we will do this operation progressively, starting from the 



root. For particle 1, we choose to put yJU\{l) on each side of U l , as in (|38|). We 
then progress along the branches of the tree; we first consider all J7 2 (l, j')'s that 
start from the root, and associate to each of them one \fUiij) that is inserted 
directly on the right side of this U 2 operator; the other y/Ui(j) is merely put at 
the end, after all the Vs. We then proceed to apply exactly the same method 
to the "second generation" of C/ 2 (j, k)'s, and again insert one ^/?7i(fc)'s directly 
on their right side, another at the end, etc.. We continue this operation until the 
whole tree is dressed with C/i's; the construction is sketched in figure ||, which 
shows the association between the additional Ui's and the initial E/2's. Applying 
this dressing procedure to each term of the right hand side of ( ^9| ) leads to an 
operator that we note Uf 1 . For instance, if I = 3, we have: 



C/ 3 «(l,2,3)= IvW) [^2(1,2) V ^M2)C/ 2 (1,3)VCM3)+ 
+ U 2 ( 1 , 3) y/vW) U 2 ( 1 , 2) X /V(2) + 



+ ^ 2 (l,2) v ^i(2)C/ 2 (2 ) 3) v ^7 1 (3) + ....j y/U^U^U^S) 

(41) 

Finally, to make sure that the operator is Hermitian, we simply define 
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it as the Hermitian part of 

[/«(1,2,...,0 = - 



CTj B (l J 2,...,0+(c?i B (l,2,... ) i)) t 



(42) 



(we will nevertheless sec that, in practice, this symmetrization has no conse- 
quence on the following calculations, which deal only with partial traces: the 
distinction between and Uf 1 is not essential here). 



3.3 Partial trace 

We now study the following partial trace with respect to particles 2, 3, I: 

T^, 3 ,...i{Ui R (l,2,3,...,0} (43) 

which is still an operator in the space of states of particle 1. Since Uf 1 is obtained 
from (|3^) by adding y/Ui's at appropriate places in each term of the sum, the 
trace contained in (^) can itself be expressed as a sum of traces of products 
of operators. To each of these terms, we will now associate another sort of tree 
diagram, which resembles that of figure |^, but where the order of the branches 
now characterize the order of operators under the trace (instead of being related 
to the numbering of the particles). 



3.3.1 Normal ordering 

The purpose of this section is to put the operators in a standard order that can 
be described by a tree-diagram corresponding to a well-defined contribution to 
the value of the partial trace that provides p\. As before, the open circle at the 
root of the tree corresponds to particle 1, but now the upper branch corresponds 
to the first sequence of operators that occur inside the trace, the first sub-branch 
to the second sequence, etc. This new tree is actually not very different from 
the initial tree: the only difference is actually the way in which the ordering of 
its branches is defined. 

Let us first consider a given term contained in the right hand side of (p9|), 
and suppose that we are interested in its partial trace - for the moment we leave 
aside the \fUis, and define the notion of "normal ordering" for the partial 
trace of any term that appear in the sum defining [7/(1,2, ...I). To reach this 
normal ordering, the first step is to locate the first operator of the product that 
contains particle 1, as well as some other particle i, [7 2 (l,z), and to move it 
to the front (the left side of the product): all operators that occurred before 
are moved to the end of the series (the right), in the same order, by using the 
property of circular permutation under the trace. The second step is to locate, 
among all operators now sitting after this operator, which is the first operator 
[7 2 (i, j) containing i, and to move it to the second position, directly to the right 
of [7 2 (l,i); this is possible since, if the intermediate operators do not contain^ 

3 They cannot contain particle i by assumption. 
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Figure 3: This figure shows a pi tree-diagram where the order of the branches 
now determines the order of operators inside a trace over variables of all parti- 
cles, except particle number 1; the first vertical branch corresponds to the first 
group of operators, the first secondary branch on the right to the next group, 
etc. Moreover, each U2 carries with it two \JTJ\ operators; in order to simplify 
the diagrams, we represent them as half-squares (rectangles). A full square 
represents a U\ operator. 



particle 1, they can be moved to the left of C/2 (1, i) and then to the end of 
the series; if they contain particle 1, they cannot contain^ particle 3, so that 
they can be moved just after 1/2(1,3). The series now begins with the product 
1/2(1, 1)1/2(1, j)- The third step is similar to the second: one locates on the right 
of 1/2(1,/) the first operator that contains particle j, 1/2(3, k), an( i moves it 
directly to the third rank, by the same method: if the intermediate operators 
contain^] neither particle 1 nor particle i, nor particle j, they are moved to the 
front and then to the end; if they contain one of them, they do not contain^ 
particle k, they are moved just after 1/2(3, k). The series now begins with the 
product U 2(1, i)U 2(1, 3)1/2(3, k). The same process continues by iteration until, 
at some point, one reaches the end of the branch of the tree, with a numbered 
particle that does not occur any other U 2 ■ 

One then proceeds to construct a new branch, and therefore to select a 
ramification point. For this purpose, one goes backwards from the end along the 
first branch, locates the first numbered particle m, in operator V 2(^x1, p), that 
occurs (at least) a second time in the list of remaining operators, in 1/2(171,11); 
this latter operator is then moved directly to the right of 1/2(171, p), and creates 
the starting point of another branch. The new branch is then extended by the 
same method as the first. When this branch is also finished, two cases may 
occur: either particle m occurs a third time (or more), so that three branches 
(or more) of the tree will originate from the same point; or particle m does not 

4 They cannot contain both particles 1 and j, since otherwise the two particles would be 
linked twice, which is contradictory with the tree structure of figure ^. 
5 They cannot contain particle j by assumption. 

6 The intermediate operators can not contain at the same time particle k and either particle 
f , or i or j: for instance, if they contained k and i, those two particles would be linked twice 
(directly and through particle j), which is contradictory with the structure of the tree-diagram. 
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occur anymore, and one continues to move backwards in the main branch to 
find another particle that occurs again in one of the non-ordered JJa's. At some 
point, all operators have been moved to their appropriate place, and the process 
stops. 

Now, to get the contribution to the trace, we have to add \/Ui , s at the 
appropriate places, but this does not change much to the reasoning that me 
have made for reaching the normal ordering of the operators: the \/TJi's at the 
end of the product do not move during the operation, while those that directly 
follow U2S move together with this operator (this is possible since they do not 
change the commutations rules). Consequently, all operators can be put into 
their normal ordering exactly in the same way, as for C/2's only. 

Obviously, the numbering of all particles, except particle 1, is totally ir- 
relevant for the value of the trace, since it defines dummy variables; what is 
important is the geometry of the tree, since it determines the value of the op- 
erator obtained after the partial trace is taken^j. Therefore, a diagram such as 
that of figure || corresponds to a well-defined contribution to the trace. Another 
remark is that the contribution of UP to the trace already provides an Hermi- 
tian operator acting on particle 1: by circular permutation under the traces, it 
is easy to see that the only change induced by a reversing of the order of the 
operators leads to a thee where the order of the branches is reversed, in other 
words to a tree that already exists in and ensures Hermiticity of the partial 
trace. This means that the Hcrmitian symmetrization of ( f42| ) is actually not 
necessary; from now on, we will therefore only consider the trace of Uf 1 . 



3.3.2 Weights of tree-diagrams 

From the results of the preceding section, we can express the partial trace of 
as a sum over all diagrams such as that of figure ^; we now wish to know 
what their weight is, or in other words, how many terms of the double sum ( |39| ) 
correspond to each diagram. The reasoning is analogous to the reasoning of § 



2.3, except that now we are dealing with order of operators, instead of particle 
numberings. Suppose that we ascribe arbitrary particle numbers to all nodes of 
this diagram - since particle 1 is always at the root of the diagram, this can 
be done in (I — 1)! different ways. To each of these numberings corresponds a 
given sequence of U2 operators, and therefore a given term of the first summa- 
tion of (|39|). Now, how many terms of the second summation then correspond 
to this given diagram? To answer this question depends on the branching (or 



ramification) factors ri, etc. defined in § 2.2 . In the list of L^'s, the only non- 
commuting operators are those which contain one common particle^. Starting 
from any order of operators, and after ordering them according to the proce- 
dure of § |3.3.1 , it is easy to see that there is a probability l/(ri)! that the first 



ramification will take place with the appropriate order of chain of operators, a 
probability l/(r!)!(r 2 )! that the right order is still obtained after two ramifica- 
tions, etc.. Finally, among all operators that are contained in the second sum 

7 But different trees do not necessarily correspond to different results. 

8 They never contain two, otherwise the diagram would not be minimally connected. 
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of (p9|), the number or terms that correspond to each tree-diagram is: 

n^)! 



(44) 



If we take into account the factor (I — 1)1 mentioned before (corresponding to 
the first summation) as well as the factor 1/(1— 1)! which is explicit in ((H), we 
obtain the following weight of the tree-diagram: 



(l-l)\ 



1 (J-l)! (Z-l)! 



so that we eventually obtain: 

rr 2l 3,...i{tfj R (l,2,3,... > 0} 



(i-iy.ik( ri )\ n^n)! 

(/-I)! 



E 



{pi diagrams} 



Tdiag.{l) 



(45) 



(46) 



where Tdiag.iX) is the operator (acting on the variables of particle 1) that is ob- 
tained by the partial trace corresponding to the diagram with branching factors 
n,r 2 , etc. 

3.4 Energy shift for p T 

The next step is to obtain the one-particle density from the previous consid- 
erations. Its expression in terms of Ursell operators || (only retaining their 
tree-reducible part) is: 



PI (1) = zUt(l) + z 2 Tr 2 {Kf (1, 2) } + Z - Tr 2>3 {U*(l, 2, 3) } 



j-i 



(47) 



+ j^ l Tr 2l .. l i{U l R (l,...,l) } + ...) + 



Replacing the reduced Ursell operators by their expression as sums of tree- 
diagrams and taking the counting factors into account (|45|) , we see that pi 
becomes: 

( □ p w \ 

2! Y 2! 



Pj(1) = ^zUx{l)ll + - xi + 



z 
21 



2! 



1! 



j-i 



(i-iy. Ui( ri y. 



V 

n n □!/□ 
u\/un\u □ 



j 



(48) 



16 



It is now possible to regroup diagrams according to the "first branching factor" 
l/(ri)!, and to write: 




z 2 h + ±zl 



(49) 



Again, we notice a self-similarity property of the expansion, which makes one- 
particle density operators appear on the right hand side of the equation: 



Pi 



(l) = y/zU : 



!eh En 
1+0+ 2j XC-HE3 + ... + ^y 




We now recognize the development of an exponential operator, and finally ob- 
tain: 



pi(l) = y/zUi(Xj exp (Tr 2 {U 2 {1, 2) Pl (2)}) y/zU^l) 



(51) 



or: 



-0A(l)=Tr 2 {U 2 (l, 2) PI (2)} 



(52) 



We therefore find that the energy shift depends linearly on the one-particle 
density operator (proportional to the density), which corresponds exactly to 
the mean-field approximation. 



3.5 Physical discussion 

The preceding calculation illustrates the physics that is behind the construction 
of a mean-field: when the test particle 1 interacts with another particle i, it 
is also possible that this particle interacts in turn with particle j, and so on. 
Moreover, either the test particle or any other particle may perfectly well inter- 
act with several others; this introduces branching in the interaction tree, either 
directly at the root or at any other place. What is not possible is to create 
"loops": any particle can interact with one or several neighbors, and these in- 
teractions can propagate further to other particles through many intermediate 
carriers, but they should never come back to the original particle. This is very 
similar to the "no re-collision" assumption which is behind the "molecular chaos 
Ansatz" of the Boltzmann transport equation. 

We note that the expression (|5^) involves the operator 1/2(1,2), instead of 
the potential itself V12 as in the usual expressions of the mean-field. In other 
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words, the mean-field involves matrix elements of an exponential containing the 
interaction potential (see definition (|3^) of the second Ursell operator), more 
precisely the difference between two exponentials (corresponding so to say to a 
change in the local Boltzmann equilibrium) , instead of not merely the matrix el- 
ements of V12 itself. Of course, if this potential can be treated to first order, this 
makes no difference, as shown by an elementary calculation. But realistic inter- 
atomic potentials can not be treated properly by first order perturbation theory 
(Born approximation), and this exponential may actually introduce an enor- 
mous difference. A well-known illustration of this fact is given by alkali atoms, 
which have a strongly attractive potential sustaining many bound states (V12 is 
negative, except in the very short range part of the potential); nevertheless, the 
phase shift at zero energy may correspond, either to a positive scattering length 
(effective repulsion), or to a negative value (effective attraction), depending on 
small details of the potential and on a very delicate balance effect between at- 
traction and repulsion^. A naive reasoning in terms of the potential itself could 
lead to the idea that, in a dilute gas of alkali atoms, the mean-field should al- 
ways be very attractive. Of course, this is known to be incorrect: the mean-field 
is actually repulsive when a is positive, and attractive only when a is negative. 
The usual way to understand this property is to replace the real potential V12 by 
a pseudo-potential, which is directly proportional to the scattering length a and 
treated to first order, a somewhat heuristic method (since the exact reason why 
using the real potential is incorrect is not so clear). Here, we clearly see that 
what appears naturally is the matrix elements of U2O-, 2); as shown in ]l9| , the 
latter can be expressed in terms of phase shifts and be shown to be proportional 
to the scattering length^ a: for instance, if the latter is positive, we directly 
get a positive value for the mean-field, without any special manipulation]^]. 

Another remark is that ( |5l| ) and (|52|) give pi as the product of exponentials, 
and not the exponential of a sum, which is not the same thing if the oper- 
ators do not commute. For translationally invariant systems, this distinction 
vanishes since all the single particle operators are diagonal in the same basis 
(the momentum basis), and therefore commute. In other cases, one should be 
careful to take into account non-commutativity of operators; for instance, the 
Hermitian operator A can not be seen exactly as a correction to the one-particle 
Hamiltonian. 

Finally, we remark that our reasoning could be generalized. Here, we have 
expressed all reducible parts of Ursell operators in terms of U2, which is its 
own irreducible part; we have left aside the irreducible part U 3 rr of U3. But 

9 Even a purely attractive potential may have a positive or negative scattering length, 
depending on the position of the last bound state with respect to the continuum. 

10 More precisely, this is true for the dominant part of the matrix element, which depends 
on the value of the collision wave functions outside of the interaction potential (asymptotic 
value of the scattering states). Another contribution arises from the wave functions inside 
the potential. Nevertheless, if the range of the potential is very small, the latter contribution 
remains negligible. 

11 Of course, other methods to prove the same result also exist; for instance, in the Green's 
function formalism, the summation of an infinite series of ladder diagrams can be used to 
construct the scattering length from the potential itself. 
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it should be possible to go further, and to use U 3 as the starting point of 
another decomposition of all Ui's for I > 4; in this way, all these Ui's could 
give a contribution to the energy shift A that is quadratic in the density. More 
details on this calculation are given in Appendix III. 



4 Quantum statistical physics: identical parti- 
cles 

The formalism for identical (Bose or Fermi) particles differs from Boltzmann 
particles by the inclusion of exchange cycles. They introduce horizontal parts 
in the Ursell operator diagrams j|, ||, |(| , which combine with the vertical lines 
associated with the Ursell operators; many more diagrams have to be taken into 
account in order to include quantum statistics. Despite this big difference, we 
will see that the formalism introduces almost the same mechanism and that the 
exponential of energy corrections also appear naturally. 



4.1 General equations 

As shown in reference ||, for identical particles, equation (|3^) should be replaced 
by: 



pi = fi 



(53) 



TH HT+W7T1 



where f\ is the one-particle density operator for an ideal gas: 



/i(l) 



zU t {l) 



1-77^(1) 



(54) 



with i] = +1 for bosons, r/ = —1 for fermions. In (53), the operators [1 + rjfi] 
arise from a summation over all possible values of cycle length I (ranging from 
I = to infinity) of the product [^zl^l)]' - the origin of the factor r\ is that, 
for fermions, every exchange cycle of length I introduces a (—1)' sign into the 
diagram which contains it. Generally speaking, in all such diagrams, the lowest 
horizontal line corresponds to the multiplication of operators acting in the space 
of particle 1, while all the other lines above correspond to traced variables of 
other particles. 

In the series giving pi was re- written in a self-consistent way, which 
resums in one single diagram an infinite number of diagrams of the initial series 
53|): 
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pi = fx + ITWIh Utw^ 



(55) 



- IT+rTFT 



For instance, the diagram of the third line of (p3|), as well as many other similar 
diagrams, are now included in the first term after fx in (55), which symbolizes 
the partial trace: 

2z 2 (1 + V h(l)) Tr 2 {U 2 (l, 2) [1 + VPl (2)}} (1 + VPl (l)) (56) 

In (^||), the second term after fx turns out to be the exchange term of (p6|), 
since it can be written: 

2z 2 (l + r,fx(l)) Tr 2 {U 2 (l,2)r,P ex .(l,2) [1 + VPi(W (1 + VPrO-)) (57) 

where P ex . is the exchange operator between particles 1 and 2. These two terms 
therefore group naturally together with the introduction of the symmetrized (S: 
rj = 1) or antisymmetrized (A: r\ — —1) form of the U 2 operator: 



tf' A (l,2) 



1 + ^.(1,2) 



C/ 2 (l,2) 



(58) 



For a dilute gas, they were interpreted in Q as the mean-field correction, which 
provides no correcti on t o the critical temperature of the gas within the MIME 
approximation (see §3.1). The explicit expression of the third term is, similarly: 



2z 2 (1 + ^(1)) Tr 2 {U 2 (l, 2) (1 + VPI (1)) (1 + npj(2)) x 
xf/ 2 (l,2) (l + W/ (2)) (1 + VPiW)} 



(59) 



and, as the preceding term, it groups with another diagram (not shown) that 
appears as its exchange diagram. 

4.2 Introducing exponentials 

Our purpose now is to introduce exponentials in order to make energy shifts 
appear explicitly. For the ideal gas, the one particle distribution is given by ( |54| ) 
so that, since rf = 1: 



Similarly: 



1+^C/ir 1 =V [zUx] 1 = V exp [/?(#! (1) -A*)] 
V [l + r)(fi)- i r 1 = ^L_=zUx 



1 + vh 



(60) 



(61) 
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Figure 4: With quantum statistics, this diagram replaces the diagram in the 
right part of figure [| note the presence of the statistical [1 + rjpi] factors. 



where the same exponential appears, now with a negative exponent. To calculate 
energy shifts, it is therefore natural to introduce quantities such as 1 + ?7(pi) -1 
or + 

To do this, we remark that the lower lines of all diagrams contained in ( p3| ) 
start and end with the same operators, so that equation (55) can be re- written 
as: 



P/(l) = /i(l) + [l + ?//i(l)] K(l) 



(62) 



1 +')PJl 



(63) 



where the operator K is an infinite sum of partial traces, containing various 
number of traced U2& as well as [1 + r)pi{i)\s in the horizontal lines: 

H-n- EElg2 jj , 1 
K = i-Uh + h- Hi+mrl -W-i 

Hr- EBflZh H 
+ J- htWpTh Ih + ... 

In other words, the diagrams corresponding to K are obtained from those cor- 
responding to pi by simply removing the two external parts of the lowest line 
(diagram amputation). The operator K is clearly Hcrmitian, as can be shown 
by using circular permutation of operators under the trace. There is an obvious 
analogy between operator K and the self-energies in the formalism of Green's 
functions. 

Now, if we multiply both sides of ( |62[ ) by (/i) _ on the left, and by (pi)^ 1 
on the right, we obtain: 



(hr 1 = (pir 1 + (1 + v k (1 + n (p!)- 1 ) 



or, if we multiply both sides by 77, add 1, and take (|6l| ) into account 
1 + V (A) -1 = [l + (^ir 1 K][l + r) (p!)- 1 ] 



(64) 



(65) 
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or again: 



Pi 
1 + VPi 



= zUi + K 



(66) 



In the absence of interaction, K vanishes and we recover (pl|); K appears there- 
fore as analogous to a correction to zU\ introduced by the interactions, in other 
words to a correction to the single particle energy. 

To make this analogy more precise, by similarity with (|60|), we define the 
energy shift operator A(l) by^ 



1 + (VPi)' 1 = V Wxil)]- 1 ' 2 [zU^l)]- 1 ' 2 
which immediately leads to: 

1 „,.s 1 



cxp{-/3A(l)} = 



1 



K{1 



(67) 



(68) 



The scries defining the operator K(l) therefore contains all the information on 
the energy shifts for the single-particle density operator; equation (|68|) is the 
equivalent of relation (pq) of Appendix II, in an operator form, when quantum 
statistics is added. 



4.3 Calculating the energy shifts 

Our purpose now is to express the energy shift as a function of the density, or of 
the single particle density operator pi, as we did in classical statistical physics. 
This raises two problems: first, we have to build an exponential from the infinite 
series that provides K; second, we have to understand how the density operator 
pi appears, instead of the factors [1 + rjpi(i)] that seem to be systematically 
present in all terms of the series. We will see that the two problems "cure each 
other" and that the result the calculation leads to an systematic and satisfactory 
grouping of the terms for a dense system. To do this, we begin with a simple 
case, where only the direct term of the mean- field is taken into account. 



4.3.1 Direct term 



We now take an approximation of expression ( |63| ) of K by retaining only the 
diagrams that contain one single U2, one single U3, ... one single Ui, ... and 
where all the Ui operators are connected to I different cycles: 

h 



K 



i+VPlI h 
+ h 



(69) 



12 We could also have used another definition of the energy shift operator, by introducing 
the exponential of a sum instead of the product of exponentials: 

It turns out that, for the present calculation, definition ( |r3^ ) is more convenient. 
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Moreover, we replace each Ursell operators Ui by its tree-reducible part. This 
leads to the equation: 



K(l)= z 2 Tr 2 {[/ 2 (l,2) [l + W/ (2)]} 

+i Tr 2 , 3 {U*(l, 2, 3) ([1 + rj Pl (2)} [1 + Wl (3)]} 

+.... (70) 

+ T ^ w Tr 2 ,..,i{U l R (l,..,l) [l + 77/>/(2)]... x [l + wiU)]} 



For the sake of simplicity, we will even replace here the operators Uf 1 by their 
non-symmetrized part Uf 1 (see equation ([l2])); nevertheless, since at the end we 
will find an Hcrmitian operator, this has no consequence on the result. Graphi- 
cally, if we represent as before the y/Ui's by vertical rectangles, Ui's by squares 
and multiplication by [1 + rjpi] by long horizontal rectangles as in figure |], ( f70| ) 
then corresponds to: 

Tnx iTTwn 
X EE5251 nx ii+ipii plnx n+ioi 

z 2 ,2\ T,^nx iTTT7pn 2! 



K(l) = Ji x6 +^xi- + -xo ) + ... (71) 

We now notice the same kind of self-similarity of the series than in the case of 
Boltzmann statistics. There are nevertheless two differences. The first is that 
the operator that is resummed in the branches is now K, or actually the sum 
zUi + K (the term zU\ is introduced by the 1 that is present at each node 
but absent as the first term in the series); the second is that this operator is 
multiplied on the right by the factor [1 + rj pi], so that we get the product: 

[zUt + K] x [1 + VPI ] (72) 

But, according to (|66|), this product is merely equal to /?/. We therefore have: 



K(l) = VCMI) [exp (Tr 2 {tf a (l,2)p,(2)}) - l] y/UJT) (73) 

(the —1 arises because the first term of the exponential series is absent from 
K). Inserting this result into ( |68| ) finally provides: 

-/JA(l) = Tra{CT2(l,2)p/(2)} (74) 

which is exactly the same result as for Boltzmann statistics. We can now check 
that the final result is Hcrmitian, so that it was indeed correct to ignore the 
Hermitian symmetrization of the operators Up'. 

4.3.2 Exchange term 

We now add exchange terms into (p9). An example of such a term is the second 
term of the right hand side of (|63|): we have already noticed that this term 
is obtained by replacing 1/2(1, 2) by U2(\,2)P ex in the direct mean-field term 
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<0 



Figure 5: Exchange diagrams for the case i = 3; the direct term from which 
these diagrams are obtained is the second term in the right hand side of (|69|). 



associated with the linear density term in K. In other words, the sum of the 
two terms corresponds to the following substitution: 



U 2 (l,2) =► U 2 (l,2) [l + r?P ex .(l,2)] 



2U S 2 A ( 



1,2) 



(75) 



where P ex .{l,2) is the exchange operator for particles 1 and 2. More gener- 
ally, the question is whether each term in the whole series of terms of (fn|), 
which eventually leads to an exponential, contains all exchange terms that are 
necessary to perform substitution (fTq) . 

In order to construct the diagrams which are exchange diagrams for the 
mean-field, it is useful to remember the origin of each Ursell diagram: it arises 
from the association of a set of Ui's (with I > 2) with permutation cycles; in 
the diagrams that we have retained so far, only one XJ\ is present, and each 
of the I particles is included in a different cycle, containing none of the other 
particles. Let us now consider two particles, i and k, that sit initially in different 
cycles; if we apply an additional exchange operator P ex (i 7 j) to the product of 
the two cycles, it is easy to see that we obtain a larger cycle with an length 
which is the sum of the two initial lengths: the two cycles merely coalesce into 
one. This immediately leads to another Ursell diagram, such as the two first 
that are shown in figure ^ (a and b) in the particular case / = 3. One can then 
repeat the operation and apply another exchange operator to the result, which 
will make two more cycles fuse together into an even larger cycle; in this way 
still another diagram is obtained, such as the last shown in figure | (c). This 
introduction of exchange may be repeated until all pairs of particles contained in 
U 2 links are exchanged, which leads to a maximum cycle of length I containing 
all of them. Altogether, 2 l ~ 1 — 1 different diagrams correspond to all possible 
ways to fuse together the various cycles; they provide all the exchange diagrams 
associated with the initial direct diagram. 

Finally, we have to replace in all these diagrams Ui by the sum of product of 
Ui operators that corresponds to as we did for the direct terms; we replace 
the "skeleton" provided by Ui by all possible trees. The result of this operation is 
the same as for the direct terms: any branching at the root introduce a product 
of operators, while branching at the other nodes introduces a product inside a 
partial trace. The situation is thus not different than before, except that two 
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Figure 6: An example of a "correlation diagram", where exchange cycles and 
Ursell operators (LVs e.g.) combine to make a "loop". This diagram provides 
a contribution which is beyond mean-field. 



sorts of links now occur, with or without exchange; the summation of the series 
then provides nothing but the exponential of a sum. We therefore obtain the 
simple result: 

-/3A(l)=Tra{tT a (l,2) [1 + r,P ex .(l, 2)] P i(2)} (76) 

which shows that it is the symmetrized form of U2 that appears naturally in 
the expression of the mean-field. 

4.3.3 Correlations 

Can we go further in this exponentiation, and try to include in the energy shift 
terms such as that corresponding to the diagram of figure ^? 

As we have seen, the exponentiation operation involves the consideration 
of operators that are the square, the cube, etc. of the lowest order operator. 
Could we somehow consider diagrams such as that of figure [| but where U 3 are 
introduced, replace U 3 by a tree diagram, and show that the result is the square 
of the initial operator? Actually, this does not seem to be possible. One reason 
is that the square of an operator such as that shown in figure |^ would involve 
particle 1 to be part of two independent cycles, a situation that never occurs 
in Ursell diagrams. For this reason, we have not been able to exponentiate 
the diagram of figure ||. One possibility is to consider this class of terms as 
smalls correction to the mean-field exponential terms, to be treated linearly (as 
was done for instance in 0). This is indeed possible above the Bose- Einstein 
transition point, but it is certainly not correct below the this point, as we 
discuss below. As a consequence, the validity of our calculations is limited to 
non-condensed Bose systems. 

4.4 Physical discussion 

The energy shift operator A(l) is reminiscent of the self-energy operator of tem- 
perature Green's function's theories, but actually it corresponds to a different 
notion. Self-energies include a notion of time, or frequency (or Matsubara dis- 
crete frequencies) , while time evolution is absent from the formalism of Ursell 
operators: only equilibrium properties are obtained, which excludes notions 
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such as time dependent response functions, collective modes and quasiparti- 
cles. In other words, only the equilibrium single particle density operator is 
calculated, which contains less information (integrated information over all fre- 
quencies) than full Green's functions. It is nevertheless interesting to see how 
the formalism manages to reconstruct the energy exponentials, with or without 
exchange cycles, and that at the end the mean-field expression simplifies to in- 
troduce the single particle density operator itself in a consistent way. It is also 
worth noticing that the relevant matrix elements to calculate this mean field are 
not those of the bare interaction potential V12 , but rather those of an operator 
1/2( 1, 2) that contains V12 in an exponential - we have already discussed in § 
|3.5| , its relation with the introduction of the scattering length as the relevant 
parameter to describe the interactions. 

Of course, one should keep in mind that the notion of mean-field is not exact: 
we have summed only limited classes of diagrams; moreover, we have replaced 
the Ui's by their reducible part £7; . In the case of bosons for instance, it is known 
[^[ ^] that correlation diagrams such as that shown in figure ^| play an essential 
role just above the Bose- Einstein transition point, and even that the theory takes 
a non-perturbative character at the transition temperature, so that no limitation 
to any finite set of diagrams (or of class of diagrams) is in principle possible. 
Below the condensation point, the situation is even more dramatic. The reason is 
that a ladder diagram with M ladders (J7 2 operators) and an intermediate state 
with an extensive population (the condensed state) is proportional to the y M ~ 2 
(where V is the volume), which makes it diverge in the thermodynamic limit 
as soon as M > 3: however small the interaction parameter is, these terms will 
always dominate the others in the thermodynamic limit. A similar situation 
occurs for the so called "bubble diagrams" . As a consequence, since we have 
ignored all these correlation diagrams, in their present form our calculations 
have a validity that is limited to non-condensed boson systems. Actually, one 
would expect that, in a condensed system, the mean-field should include two 
different parts, due to excitations and to the condensate respectively, but for 
the moment we have not explored this question. 



5 Conclusion 

At the end of their 1938 article, Kahn and Uhlenbeck remark that "[they] 
have not been able to generalize this physical interpretation of the (3i [(irre- 
ducible cluster integrals)] to the quantum theory." In the present article, we 
have achieved this goal: we have provided a consistent derivation of the expo- 
nentials that introduce the energy shifts A, from which the equation of state 
can in turn be derived as in classical statistical physics. Actually, in the article, 
we emphasize how the linear density term can be derived, but Appendix III 
discusses briefly how higher order density terms could also be included. For a 
classical system, our method is different, and in a sense simpler, than that of 
Hansen and McDonald jig ]; it does not require any reasoning in the complex 
plane and we have shown how the exponentials of pi appear rather naturally in 
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the calculations. The operator that plays the basic role in all our calculations 
is 1/2(1, 2); its matrix elements depend on the asymptotic (long distance) prop- 
erties of the two body scattering wave functions. But, in fact, 1/2(1, 2) contains 
more information than only asymptotic wave functions and phase shifts: it also 
contains information about short range effects of the potential (atoms in the 
middle of a collision) as well as about bound states; it would be interesting to 
explore their consequences on the properties of the mean-field. 



The Laboratoire Kastler Brossel (LKB) is Unite Mixte de Recherche du 
CNRS (UMR 8552) et de PUniversite Pierre et Marie Curie (Paris). 



APPENDIX I 

In this appendix, we explain why classical tree-diagrams consisting of I par- 
ticles all have the same numerical value, namely (/3i) /_1 .Let us consider the 
following integral corresponding to a I particles tree-diagram: 



T d ia g . = y J d 3 r 1 d 3 r 2 ...d 3 n fu-fij -fkl 



(77) 



Each Mayer function appearing in it depends on the relative distance between 
two particles: 

f ij = f(\r i -r j \) = f(r ij ) (78) 

Let Xfc be the relative position vector between particle k and the preceding 
particle in the branch starting from the root. There is Z — 1 such vectors (k = 2 
to I). The important thing to notice is that each Mayer function depends on 
only one Xfe, and that two Mayer functions in the product necessarily depend 
on different X&. Let 

R= y(ri+r 2 + ... + r ( ) (79) 

be the center of "mass" position vector. Then making a change of variables 
in the integral from (n, r2,..., r/) to (R, X2,..., x;) (the Jacobian is one), we 
obtain: 

T diag . = y J d 3 R d 3 x 2 ...d 3 Xl f(x 2 ) ...f(xi) (80) 

We can then integrate over the variable R to get rid of the volume factor. The 
I — 1 integrals that are left separate so that: 

Tdiag. = J d 3 x 2 f(x 2 )...J d 3 Xl f(xi) = d 3 x f{x)j (81) 
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But the integral of the Mayer function is precisely what we called the first 
irreducible cluster integral f3\ (see equation (Il5])), so that we finally get: 



T diag . = Wi) 1 - 1 (82) 



APPENDIX II 

In this appendix, we emphasize the similarity between the equation that gives 
the energy shift as a function of the density and the elimination of z between 
(||) and (^) that provides the equation of state. For simplicity, we assume that 
no external potential acts on the particles, so that translational invariance is 
satisfied. We call pk the diagonal elements of pj : 

p k =< k | pi | k > (83) 

and iti (k) those of U±: 

ui(k) =<k|C7i |k>=e - " e » 

where is the single particle kinetic energy. As for the diagonal elements of 
XJi in the momentum representation, we note them: 

w 2 (ki,k 2 ) =< ki,k 2 I U 2 I ki,k 2 > 

(for the moment we do not make the MIME approximation), u 3 (ki,k 2 ,k 3 ) for 
those of /7 3 , etc. Equation p|) then provides: 



(84) 



p k =tti(k) z + z 2 J2 k ' U2(ki,k 2 )ui(k') + 
+7Ek' Ek" u 3(ki,k 2 ,k 3 ) u 1 (k')u 1 (k") 



Now, if we introduce the energy shift A(k) by the relation: 

pu = zu 1 (k) xe-W (85) 

we obtain: 

e -/3A (k ) =i + 2! ^; k/U2 (k 1) k 2 )ui(k')+ 

+4 E k ' E k " « 3 (ki,k 2 ,k 3 ) Ul (k') Ul (k") + ... 



(86) 



In this relation, we express the energy shift A(k) as a function of a series in 
power of the populations z ui(k) of the ideal gas; but if we write: 

e -/3A(k) =l + Ek , U2( k 1) k 2 )e^( k ')p k , 

,lv, , , s /3[A(k')+A(k")l 

+ lEk' Ek" w 3 (ki,k 2 ,k 3 ) e L J P k 'Pk"+-- 
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we obtain another expression where only the actual populations p k ' play a role, 
and from which the variable z has now disappeared. Finally, it is sufficient to 
take the logarithm of the right hand side of this series to obtain the expression 
of the energy shift A(k) as a function of the populations as well as all the other 
energy shifts A(k ). 

In the MIME approximation, the interactions constants U2, W3, etc. factorize 
out of the sums, the energy shift becomes a constant A that is independent of 
k, and one merely obtains: 

e-^ = l + u 2 eJ jA N+^e 2 ^N 2 + ... (88) 

where N = (N) is the mean total number of particles (the trace of pi). This 
equation provides an implicit equation between the energy shift A and the 
density p = N/V (since 1*2 is proportional to the inverse volume, 113 to the 
square of this inverse volume, etc.) From this result, one can expand the energy 
shift in powers of the density: 

A = u 2 N + u' 3 N 2 + U4 iV 3 + ... (89) 

where the new coefficients it 3 , u 4 , ... can be calculated step by step by inserting 
this relation into ([S8|); here again, the first term in the right hand side of ( |S9"| ) 
corresponds to the mean-field, and the following terms to density corrections. We 
note the similarity between this result and the usual elimination of z between the 
two z-expansions of the pressure and the density, which provides the equation 
of state, where only the actual density of the system appears. We also remark 
that, in general, the energy shift has a contribution that is linear in density, 
which can be called mean-field, but also higher order density terms (density 
corrections to the mean- field) . 



APPENDIX III 

We have defined in section the reducible part U l (1, 2, I) of any Ursell 
operator {7/(1,2, ...V) in terms of "trees" made of chains of U2 operators. For 
instance, the Ursell operator of order I — 3 then becomes the sum of a tree- 
reducible part and of an irreducible part: 

1/3(1,2,3) =Z7 3 fl (l, 2, 3) + LT 3 >r -(l,2,3) (90) 

In this appendix, we sketch how the notion of reducibility could be generalized: 
the irreducible part U 3 could be used, exactly as U2, as a starting point to 
define a "t/3-reducible" part of all operators Ui with I > 4 in terms of products 
of U 2 s and U 3 rr ''s. 

The basic idea is to build trees, not only with lines that represent C/2's, but 
also with triangles that represent U 3 's. Figure |7| shows an example of such a 
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Figure 7: An example of a diagram including irreducible parts t/ 3 of the three 
particle Ursell operator t/3, symbolized by triangles, in addition to straight lines 
representing the U^s. Such diagrams contribute to a tern in the energy shift 
which is quadratic in pi. 



tree: at each branching point, one may now connect either lines, or one corner 
of a triangle; in the latter case, the two other corners can be used as starting 
points in order to extend the branch further, with any product of t/2's and 

ITT Iw 

U 3 's. One difference is that, because each U 3 , once added, introduces two 
new possible branching points (instead of a single one for at/2), the question 
arises as to which one is represented first in the clockwise order; it is easy to 
see that this introduces a 1/2 factor per U 3 r contained in the diagram into 
the corresponding weight. Otherwise, not much is changed, provided of course 
all branching factors (including those resulting from t/ 3 's) are included in the 
weight; in other words, in formula (p8f), a factor (1/2)™ 3 is added (where n 3 is 
the number of t/ 3 rr ''s), and the ramification factors n now include the effect 

of t/ 3 's. This being done, the property of self similarity appears again in the 
series, and the essence of the reasoning still holds with relatively minor changes. 
Finally, one is led to the introduction of an exponential: 

e -^ A « (91) 

where A(l) would now be given by the sum: 

-0A(1) = Tr 2 {t7 2 (l, 2)^(2)} + 7> 2 ,3 [^Ui rr (l, 2, 3^(2)^(3)} (92) 

In other words, the energy shift would no longer be proportional to the single 
particle density operator, but would contain a term that is quadratic in pj. 

Similarly, one can expect that the irreducible part t/ 4 rr of t/4 would lead to a 
contribution to the energy shift that is cubic in pi, and so on. We have not yet 
performed these calculations. 
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